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^ . Abstract 

We analyse the homogeneous instabihties in a nematic hquid crystal subjected 

^ I to plane steady Couette or Poiseuille flow in the case when the director is pre- 

O 

aligned perpendicular to the flow plane taking into account weak anchoring at 
^ ^ the conflning surfaces. The critical shear rate decreases for decreasing anchor- 

fS| ' ing strength and goes to zero in the limit of torque-free boundary conditions 

For Poiseuille flow one has two types of instability depending on the values 

T 1 

I of the azimuthal (Wa) and polar (Wp) surface anchoring strengths. The crit- 

I ical line in (VFa^W^p) space which separates the two instabilities regimes is 

, obtained. 

o 
o 

^ I. INTRODUCTION 

X. 

■ During the last decade the study of surface anchoring in nematic liquid crystals (NLCs) 

for different types of confining substrates has attracted much attention. In particular, a 
change of anchoring strength strongly influences the orientational behaviour and dynamic 
response of NLCs under external electric and magnetic fields. This changes the switching 
times, which play an important role for applications (see, e.g., The hydrodynamic flow 
is a crucial ingredient for the dynamic response and switching characteristics of liquid crystal 
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devices. 

The anchoring characteristics can also be studied in orientational phenomena induced by 
hydrodynamic flow. In fact, recently the surface orientational transition caused by oscillatory 
shear flow was found [Q] and the influence of weak anchoring on the linear response of 
the NLC to oscillatory flows was studied Up to now the study of orientational bulk 
instabilities in NLCs under hydrodynamic flow have been restricted to the case of strong 
surface anchoring (fixed director orientation at the confining plates) 0-^]. 

In this paper the influence of surface anchoring on the homogeneous instabilities in NLCs 
subjected to steady flow of Couette and Poiseuille type is studied theoretically in the case 
when the director at the bounding plates is oriented perpendicular to the flow plane. This 
is the simplest geometry because the initial state with the director oriented everywhere 
perpendicular to the flow plane is, by symmetry, a solution of the nematodynamic equations 
for any shear rate. The state can change only via a symmetry-breaking instability. The type 
of instability strongly depends on the sign of the Leslie viscosity coefficient a^. The case of 
negative 03 {flow- aligning materials) and strong anchoring of the director was investigated 
theoretically by Leslie 0] and Dubois- Violette and Manneville PJ^-[I(]|. They have shown 
that in the absence of external fields the first instability is a homogeneous one in both, steady 
Couette [0,11 and Poiseuille p fiows . The theoretical results were found to be in good 



agreement with experiments of Pieranski and Guyon and coworkers ||Tl|--|13| . In Couette fiow 
sufficiently strong additional magnetic field applied parallel to the initial director orientation 
was found to change the type of instability into a spatially periodic one where rolls develop 



rT|,|T3|, in contrast to the case of Poiseuille fiow, where a magnetic field does not change the 



type of instability [|T^ . Well above the threshold of the homogeneous instability in Poiseuille 



fiow rolls were observed to develop in a secondary instability In the case of positive 
{non- flow- aligning materials), according to the mechanism of Pieranski and Guyon, one has 
no homogeneous instability and only rolls are expected ||TT|. Since for these materials there 
is also a tumbling motion [|r^ , the orientational behaviour of NLCs can be quite complex. 
Here we focus on flow- aligning nematics. Starting from the standard set of nematody- 



namic equations in the Leslie-Ericksen formulation the basic equations describing the 
homogeneous instabihties are presented (sec. |I|) taking into account arbitrary surface an- 
choring strength. Full semi-analytical solutions together with approximate expressions for 
the critical shear rate at which the homogeneous instabilities develop are obtained for steady 
Couette (sec. fT|) and Poiseuille (sec. [TVD flow. In sec. |V] we discuss the results. 

II. BASIC EQUATIONS 

We consider a NLC layer of thickness d sandwiched between two parallel infinite plates. 
The origin of a Cartesian coordinate system is placed at the centre of the layer with the z 
axis perpendicular to the bounding plates. Steady Couette flow is generated by one plate 
(at z = d/2) moving with constant velocity along the x direction and the other plate 
(at z = —d/2) fixed. Steady Poiseuille fiow is obtained when a constant pressure gradient 
AP/ Ax is applied along the x axis. The confining plates provide a director orientation along 
the y axis. The basic state is given by the stationary homogeneous solution of the standard 
set of nematodynamic equations (Leshe-Ericksen equations |^) 

n° = (0,1,0), v° = (t;°,0,0), (1) 

where 

= 1/0(1/2 + z/d) - for Couette fiow , 

= -1^^(1/4 -^Vrf^) - for Poiseuille fiow, 
and ?73 = 04/2 with 04 the corresponding Leslie viscosity. 

In order to investigate the stability of solution (|I]) we linearise the nematodynamic equa- 
tions 0] with respect to perturbations that are homogeneous in the plane of the layer 

n = n° + (ra^,ny,n^) , y = \° + {vx,Vy,Vz) , (3) 

where rij, Vi [i = x, y, z) are functions of z. We are looking for the existence of a stationary 
solution of the linearised nematodynamic equations, which signals the onset of a stationary 
(i.e. nonoscillatory) instability. The linearised equations are 
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0!2Snz - K22nx^,, = , (4) 

asSn^ - Kiin^ ^^ + a^Vy ^ = , 

where S = ^ is the shear rate, r]i = {a^ + + aQ)/2 and Kn, K22 are the elasticity 
coefficients for "splay" and "twist" deformations, respectively. The notation /^^ = df/dz is 
used throughout. 

The boundary conditions for the y-component of the velocity perturbation are 

Vy{z = ±d/2) = . (5) 

The surface anchoring of the director is described in terms of a surface energy per unit area 
Fs which has a minimum when the director at the surface is oriented along the easy axis 
(parallel to the y axis in our case). It is convenient to characterise the surface anchoring by a 
"polar" anchoring strength Wp pertaining to out-of-substrate-plane director deviations and 
an "azimuthal" anchoring strength Wa related to distortions within the substrate plane. 
A phenomenological expression for the surface energy Fg can be written in terms of an 
expansion with respect to (n — n°). For small director deviations from the easy axis one 
obtains for the surface energy 

Fs^^Wanl + ^Wpnl, Wa>0,Wp>0. (6) 

The boundary conditions for the director perturbations can be obtained from the surface 
torques balance equation 

OF OF 
±i^22n.,. + ^ = 0, ±Xnn,,, + ^ = 0, forz = ±d/2. (7) 

Introducing the dimensionless quantities 



5 = ^, S^PraS, ra^^=^ 



■2i, 



with 



(8) 



z — "^z 



2^03^1 

"2 Kii b ?73 



Eqs-d^) can be rewritten in the form 



:i~b){SN, 



. 



For the shear rate S one has 



The boundary conditions (|^), (0) reduce to 



for Couette flow , 
for Poiseuille flow 



± iV, ^ + WaN, = , ±Ar,, + WpN, = 0, K, = 0,atz = ± 1/2, 



where 



Wad/K22 



(9) 



(10) 



(11) 
(12) 



(13) 



(14) 



In the limit of strong anchoring {wa, Wp oo) one has = Nz = aX z = ±1/2 whereas 
for torque-free boundary conditions {wa, Wp — > 0) N^^^ = N^^^ = at the boundaries. From 



(|l^) one can see that by changing the thickness d, the dimensionless anchoring strengths Wa 
and Wp can be varied, the ratio Wa/wp remaining fixed. 

Solving the system of linear ordinary differential equations ( p!OD with boundary conditions 
(p!3|) one can obtain the critical value of the shear rate a^, at which the initial state loses 
stability, and determine the influence of anchoring strengths Wa and Wp on a^. 



III. COUETTE FLOW 



In this case the shear rate 5* is independent of z (0), so that the system (|T0|) can be 
solved using standard theory of ordinary differential equations with constant coefficients. 
Further, the {z —z} symmetry of Eqs. ([lOD allows for two possible types of solutions: 



l:{Vy - even; , - odd} - "odd" solution , 
II: {K^ - odd; A^^' , Nz - even} - "even" solution . 

We will always classify the solutions by the z-symmetry of the x-component of the director 
perturbation N^. For the odd solution one obtains 

Vy = Ci(l — b)acosh{az) — C2(l — b)acos{az) + C3 , 

Nx = Ci sinh(az) + C2 sin(az) , (15) 

Nz = —Ci sinh(a2;) + C2 sin(az) . 



Taking into account the boundary conditions ([13|) the solvability condition for the Ci 
("boundary determinant" equated to zero) gives the expression for the critical shear rate 

2waWp sinh(ac/2) sin(ac/2) + 

o-duia + Wp)[cosh(ac/2) sin(ac/2) + cos(ac/2) sinh(ac/2)]+ (16) 
2al cosh(ac/2) cos(ac/2) = . 

For the even solution one obtains 

Vy = Ci(l — fe)asinh(az) + C2(l — b)asm{az) — C-^a^hz , 

Nx = Ci cosh(az) + C2 cos{az) + C3 , (17) 

Nz = —Ci cosh(az) + C2 cos(a2;) . 

The boundary determinant condition now gives 
a|^6sinh(ac/2) sin(ac/2) + 

^alb{wa + Wp)[cosh.{ac/2) sin(ac/2) — cos(ac/2) sinh(ac/2)] — 
{ac6cosh(ac/2) cos(ac/2) + 
(1 — 6)[sinh(ac/2) cos(ac/2) + sin(ac/2) cosh{a c / 2)]} WaWp = . 

Note, that the expressions (p^ ) and ( [T8| ) are both symmetric under exchange of Wa and 



(18) 



Wp. This results from the fact that A^^^ and A^^ have the same symmetry. We found that 
for MBBA material parameters at 25° [|l^ {b = 0.58) the critical shear rate for the odd 
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mode is systematically higher than for the even mode, so there is no transition between 
them under variation of the surface anchoring. For strong anchoring {wa = Wp = oo) the 
expression for the critical shear rate of the even solution ( [T8| ) recovers the result obtained 
by Leslie 0. Weak surface anchoring reduces the critical shear rate compared to the case 
of strong boundary conditions. In the limit of one of the surface anchoring strengths going 
to zero {wa ^ or Wp 0) one has — * 0. In figure [l| the critical shear rate for the 
even solution calculated from Eq.(|l^) using the material parameters of MBBA (6 = 0.58) is 
shown as a function of l/wa, l/wp. 

[Insert figure ^ here] 

In order to obtain an easy-to-use expression for the critical shear rate of the relevant 
even mode as a function of the surface anchoring strengths one can use the single-mode 



approximation in the spirit of a Galerkin expansion |18|. We choose 



Vy = Ci sm{27iz) , 

= C2Kcos(7r2;) + 7r] , (19) 

Nz = C3[WpCOs{TTz) + n] , 



which satisfy boundary conditions (plSf). Substituting this ansatz into (|T^) where 5* 
and solving the algebraic system for Ci, {i = 1, 2, 3} obtained after projecting the first of 
Eqs.([TO|) onto the first mode of Eqs. ([T9|) , the second equation onto the second mode etc., 
one gets 

flc = STC^^WaWp{4: + Wa){4: + Wp)/{piP2) , 

Pl = WaWp + A{Wa + Wp) + 277^ , (20) 

P2 = lQwaWp{l - 6) + dn'^bpi . 

The single-mode approximation ( pO]) gives almost the same values for as (|T8D (relative 
error is smaller than 1%), so that they could not be differentiated in figure |l|. 
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IV. POISEUILLE FLOW 



In the case of Poiseuille flow the shear rate S is given by Eq.(12) and the symmetry 
properties of system (|1^) give us the following two types of solutions: 

l:{Vy,N^ - odd;iV^ - even} - "odd" solution , 
II: {Vy, - even; A^^ - odd} - "even" solution . 

The odd solution corresponds to the splay-mode and the even solution to the twist-mode 
in the notation of Dubois- Violette and Manneville Integration of the first equation in 
system (JT^) gives 

Vy^^ = K-{l~b)a^zN,, (21) 

with K = ioT the even solution. For the odd solution the integration constant K must be 
nonzero (except for 6=1, see below). After eliminating Vy ^ from the third equation in (p!OD 
one arrives at 

a^zN, -N^,^ = , 

(22) 

a^N, - iV,,,, = K . 

Following [0 we perform the transformation and introduce new variables lA = 

+ N^,V = N^- N^, leading to 

Uzz -ZU = -Ka-^l^ , ^ ^ 

(23) 

The general solution of (p3| ) can be expressed in terms of Airy functions Ai{Z), Bi{Z), 
Gi{Z) ^. 

Let us first consider the case of the even solution {K = 0). Then one has 

= Ci[Ai{Z) + A%{~Z)\ + C2\Bi{Z) + Bi{-Z)\ , ^ ^ 

(24) 

= Ci[Ai{Z) - A%{-Z)\ + C2\B%{Z) - B%{-Z)\ , 
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where Ci, C2 are constants to be determined from the boundary conditions (0). Substitu- 
tion of ( |2^ into (p!3| ) gives the criterion for the threshold shear rate 



WaBi^ + aV^Bi\ WpBi + aV^Bi ^ 



(25) 



Here 



= M ± M {-\all') , B^^ = B^ (^a^/^) ± B^ [-\all') (26) 

and similar for Ai\ and -Bi^. The limit of strong anchoring [wa 00, Wp 00) gives the 
result of Dubois- Violette and Manneville 

Bi{\aT) Bi{-\aT) ' 

leading to a1 = 102.59. The fact that and have different z-symmetry leads in Eq.(p5D 
to an asymmetry in the dependence of the critical shear rate on Wa, Wp, in contrast to the 
case of Couette flow. The critical shear rate [Eq.(p5D] retains a finite value when the polar 
anchoring strength Wp (which mainly controls A^^; perturbations) vanishes, whereas a"^ —>■ 
if the azimuthal anchoring strength Wa 0. 

For the odd solution with nonzero K the solution of (|2^) has the following form 



(27) 



= kK{CiAi'{Z) + C2Bi-{Z) + Gt-{Z)} , 
N, = kK{CiAt+{Z) + C2Bi+{Z) + Gt+{Z)} , 

where 

k = Ai^{Z) = Ai{Z) ± Ai{-Z) etc., for Bi^{Z) and Gi^{Z) . (28) 

The coefficients Ci in ( p7|) are determined from the boundary conditions (0). Integrating 
Eq. (^) and taking into account the boundary conditions Vy{z = ±1/2) = one obtains 
the expression for the critical shear rate of the odd mode 



/•1/2 

K - (1- b)al / zN:^{z; al, Wa, Wp)dz = . (29) 
J-1/2 
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Since is proportional to K, this undetermined integration constant falls out from (p9D. 
From thermodynamical conditions P,pO| the parameter h has to be positive, but is otherwise 
not restricted. The point 6=1 requires special consideration. In this case, from Eq.(pTD 
and boundary conditions for V^, it follows that K = {] and = 0, so that one should 
solve Eqs.(P^ keeping the symmetry of A^^; and corresponding to the odd solution. In 
other words, if 6 = 1 then one has the homogeneous instability of odd type with zero velocity 
perturbation, like in the case of Freedericksz transition. Moreover, in the case 6=1 Eqs.(p2|) 
became invariant with respect to change {A^^. -v^ A^^}, so that the critical shear rates for even 
and odd solutions are the same up to transposition {wa -v^ Wp}. 

The instability of the odd mode is mainly controlled by the polar anchoring strength Wp. 
In the limit of zero azimuthal anchoring strength — one has a finite critical shear rate, 
whereas — for the polar anchoring strength — > 0. 

The critical shear rate a1 for the even solution [Eq.(]25|)] and odd solution [Eq.(p9D] have 
been calculated for the material parameters of MBBA for various values of Wa-, Wp (figure 
Depending on the azimuthal and polar surface anchoring strengths one can have different 
z-symmetry of the first unstable mode. The critical line in (w^, Wp) space corresponds to 
the crossover between the two types of unstable solutions. 

[Insert figure here] 

Since the expressions ( ^51) and (^) for the critical shear rates of the two possible unstable 
modes are quite complicated we use a single-mode Galerkin approximation in order to obtain 
easy-to-use formulas. 

For the even mode we can use Eqs. (p2D with K = {]. Assuming 

= Ci[w a cos{n z) + it] , 

(30) 

A^^ = C2[wpSm{2Trz) -\- 2tt sm{nz)] , 
after substituting (|30| ) into ( p2D and projection one obtains 

al = QVSn^^WaiA + wj(37r2 + 2Qwp + 3wl)/pl , 

(31) 

Pl = 16 WaWp -\- 9lT'^{Wa + Wp) -\- 7 27!-"^ . 
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For the odd mode we choose 



(32) 

Nz = C2[WpCOs{7Tz) + n] . 



Substituting this ansatz into ( ^21) one finds the coefficients Ci and C2. Then from (^) one 
obtains 



al = QV3TTywpiA + Wp)i3TT^ + 20wa + 3wl)/ipiP2) 



(33) 

P2 = bpi + (1 - 6)(97r2 - 144 - 2wa)wp . 

The approximate expressions (^) and (^) give systematically higher values for the critical 
shear rate compared to (j2^) and (^), about 10% for both modes. Equating ( PT| ) and (^) one 



can obtain an approximate expression for the critical line in the {wa, Wp) plane corresponding 
to the crossover between the critical even and odd modes (see figure |^). We now write 
Pa = ^/wayPp = 1/Wp. For strong azimuthal anchoring {Pa = 0) the transition from even 
to odd mode takes place at some critical value of polar anchoring strength: Pp = PpQ, where 
Ppo is the solution of the algebraic equation 

97T%P'^, + 2n\39b-19)P^^,+ ^^^^ 
[-9n\l -b) + 1206 - ^]ppo - 18(1 -b)=0. 

For MBBA material parameters one finds Ppo = 0.307. Assuming the deviation 0a, Pp) from 
(0,/5po) to be small one has up to first order in Pa-, Pp for the critical line 

Cl{PpQ + Pp)+C2pa = ^. (35) 

where 

ci = 27bTx^PlQ + 47r2(396 - 19)/3po - 97r2(l -b) + 1206 - ^ , 

C2 = 10867r4/33o + 7r2(967r2 + 9366- 680)/52o- 36(37r2 + 40)(l-6)/5po+ (36) 
+ 2166- ^ . 

The results for the critical line for the crossover between even and odd critical modes obtained 
from the numerical solutions of Eqs. (|25|) , ( p9|) and approximation (^) are shown in figure |^. 
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[Insert figure here] 

Note, that (3a = '^/wa, (3p = 1/wp are inversely proportional to the thickness d of the NLC 
layer. Therefore, varying the cell thickness, one can cross the critical line separating the two 
regimes. For that purpose it is necessary to have the ratio Pp/ Pa = KuWa/ {K22Wp) larger 
than the slope of the critical line (~ — C2/C1). Using the material parameters of MBBA one 
obtains Wa/Wp > -C2K22/ {ciKu) = 2.28 ( Wa/Wp > 1.78 from numerical results). For 
small values of Wa and Wp, the crossover is given by Wa = Wp/b. As was said above, if 6 = 1, 
Eqs.([2^) are invariant with respect to {A^^: N^}. This means that the crossover line for 
6 = 1 is exactly defined by Wa = Wp. The same result follows from approximation formulas 
(M) and (M). 



V. CONCLUSIONS AND DISCUSSION 

It was found that changes of the anchoring strengths can cause a crossover between two 
types of homogeneous instabilities induced by steady Poiseuille flow, in contrast to the case 
of steady Couette flow, where the first unstable mode is always the even one. Semi-analytical 
expressions for the critical shear rates for both Poiseuille and Couette flow are presented, 
together with simple approximate formulas of good accuracy. 

The effect of the reduction of the critical shear rate for homogeneous instabilities under 
weak surface anchoring can in principle be used for the determination of the polar (Wp) and 
azimuthal (Wa) anchoring strengths. For that purpose one may measure the critical shear 
rates for two cells with different thicknesses and use the full expressions (or approximate 
formulas) to calculate both Wp and Wa- Alternatively, one can use an additional electric 
field applied across the NLC cell. Depending on the sign of the dielectric anisotropy of the 
NLC the electric field will stabilise or destabilise Uz fluctuations and affect the critical shear 
rate differently for different values of the polar and azimuthal anchoring strengths. This 
study is in progress. The advantage of using measurements of the critical shear rate for the 
determination of polar and azimuthal anchoring strengths compared to the other methods 

12 



(Preedericksz transition, orientational transition with hybrid orientation, small-angle light 
scattering) is that here one does not need to modify (or rebuild) the experimental setup in 
order to obtain Wp and Wa at the same time. 

It would be particularly interesting to investigate experimentally the orientational be- 
haviour of NLC under steady Poiseuille flow for a cell with Wp, Wa close to the crossover 
line separating critical modes of different symmetry. 
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FIGURES 



FIG. 1. Critical shear rate of Couette flow vs. anchoring strengths, b = 0.58 (MBBA). 



FIG. 2. Critical shear rate of Poiseuille flow vs. anchoring strengths, b = 0.58 (MBBA). 



FIG. 3. Critical line separating odd and even regimes in Poiseuille flow: numerical solution 
(solid) and approximation (dotted). 6 = 0.58 (MBBA). 
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